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Abstract 

In this paper the systematic method of dealing with the arbitrary decorations of 
quasicrystals is presented. The method is founded on the average unit cell formalism 
and operates in the physical space only, where each decorating atom manifests itself 
just by an additional component of the displacement density function in the average 
unit cell. Such approach allows us to use almost all classical crystallography algo- 
rithms for structure refining based on experimental data and may meaningly decrease 
the number of parameters which have to be fit. Further help for such analysis may 
be the use of proposed recently average Patterson function, here applied to decorated 
sets. As an example we present a description of a class of decorated quasicrystals 
based on Sturmian sequence of two interatomic spacings: we calculate explicitly 
structure factor, the shape of average Patterson function and give an algorithm for 
pattern analysis. 

PACS numbers: 61.44.Br, 61.43.-j, 61.10.Dp 

Keywords: quasicrystals, disordered solids, theories of diffraction and scattering 

1 Introduction 

Discovery of quasicrystals opened a completely new world for crystallographers but in the 
same moment left them almost unarmed - all sophisticated and powerful methods devel- 
oped to deal with strictly periodic structures appeared to be limited and insufficient for the 
description of this new long range atomic orderings. Elegant mathematical constructions 
gave us deep insight into these structures, but they scarcely very handy when one has to 
deal with real data obtained from scattering experiments. We focus here on problem of 
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decorations. Obviously there is no theoretical difficulty to analyze arbitrary decoration us- 
ing standard projection techniques, nevertheless the most common usage of atomic surfaces 
has to deal with their extremely complicated shapes and numerous fittable parameters. 

In this paper we propose a technique, which by operating solely in the physical space 
opens the possibility of reducing the number of parameters which have to be fit. Our 
method originates from the concept of reference lattice suggested in pQ. Each of the 
decorating atoms appears as an additional component in so called average unit cell, which 
has very similar interpretation to that of the unit cell for periodic crystals; this allows for the 
utilization of the classical crystallography refining techniques in case of quasicrystals. As 
a useful tool we propose generalization of average Patterson function approach (proposed 
recently in [2]) to the case of decorated structures. Description of quasicrystals using 
multidimensional Patterson analysis has been already developed in [3], but our approach 
is rooted in different formalism. Our work may be understood as generalization of 0]. 

The paper is organized as follows. First a short survey through the average unit cell 
formalism is given, which is subsequently used for systematic description of decorations. 
Using the new tool we give an algorithm of a structure refining for a class of ID quasicrystals 
founded on Sturmian sequence of two interatomic spacings. The work ends with specific 
example. 

2 The formalism 

The concept of the reference lattice has been already proposed in |IJ, here we present only 
only the key ideas and features which we are going to utilize in the paper. 

Suppose that A is an orthogonal projection of Delone set A' C MJ 1 , where n is arbitrary, 
onto line I: 



so x n G A is a sequence of collinear points lying on /. A may constitute Delone set itself or 
even periodic crystal. 

If there exist n-dimensional Fourier transform of set A' there exists also ID Fourier 
transform of a set A, being essentially a cut through a /c-space of A' along direction I. The 
intensity measured along this direction is: 



where f n is a scattering power of the given atom in A'. 
2.1 Structure factor 

In the average unit cell formalism for one type of atom in the system (with scattering 
power /) normalized intensity can be written as 




(1) 




(2) 



n=l 



I{k) = \F{k)\\ 



(3) 
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where the structure factor 



F (k) = f lim — V exp(ikx n ) = / Km — Vexp(i/cu n ) = / / P x (u) exp 



iku)du. 



(4) 



We call the series u n the displacements sequence of x n (induced by the reference lattice 
with period A = 2ir/k): 



u, 



X ., , 



A 



A. 



where ||-|| is the nearest integer function: if m G Z and m < x < m + 1 then: 



\x\ 



m if x G [m, m + 1/2), 
m+1 if x G [m + 1/2, m + 1). 



Any series such that 



U, 



u, 



A 



A. 



(5) 



(6) 



(7) 



will be called an unreduced displacements sequence (of x n ) if its elements belong to any 
finite interval. 

Statistical distribution P\(u) of distances u n will be called called displacement density 
function (DDF). All considered DDFs are normalized to 1 and have value outside interval 
[— A/2, A/2]. Statistical distribution P\{u) of u' n is called extended displacement density 
function (EDDF) for P\{u)\ we have: 



Px{u) = ^P A (n + mA). 



(8) 



P\(u) function may be supported on the finite interval different from [— A/2, A/2], but it is 
still normalized. Note that we can use P\(u) instead of corresponding P\(u) in all Fourier 
integrals like (@J). In general for each k we should construct different P\(u). 
There is a handy theorem we will utilize in future, it comes from j^J. 



Theorem 1. Let x n = a n + f3 n be a sum of two real series. If d n is a displacements 

(9) 



sequence of a n induced by a reference lattice, then 

u' n = d n + p n 

is an unreduced displacements sequence of x n induced by the same lattice. 
Proof. Let A be the period of the reference lattice. We have to show that 





\ = u' n - 
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A 




A 



A. 



(10) 
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Note that for any real number x, y, \\x - 
We can write the right hand side of (fTU|) as 



\x\ 



\y\\, since \\y\\ is an integer. 



(3 n -\- 0L n 



A 



A 



Pn + a„ - \\a n /X\\ A 



A 



A 



A. 



which is the left hand side. 



□ 



The structure factor for mk, m G Z can be computed from the reference lattice for k] 
it is the m-th Fourier mode of the distribution P\(u). Thus, a single average unit cell is 
sufficient to analyze structures whose scattering occurs at multiples of a fixed scattering 
vector k . This situation includes, but is not limited to, the case where the original point 
pattern was periodic with period 2n/k . 

For modulated structures (including quasicrystals), there are usually two periods, ao 
and 60, which may be incommensurate. Using two reference lattices, the first one having 
periodicity ao an d the second having periodicity bo, the structure factor for the sum of two 
scattering vectors k = 2n/a and q = 2ir/b can be expressed by: 

1 N 1 N 

F{k + go) = hm — exp(i(A; + q )x n ) = lim — exp(i(k u n + q v n )) 



71=1 

a /2 f b /2 



n=l 



P ao ,b ( u i v ) exp(i(k u + q v))dvdu, 



(12) 



a /2 J -bo/2 



where u and v are the shortest distances of the atomic position from the appropriate points 
of two reference lattices and P ao ,b {u, v) is the corresponding probability distribution, which 
thus describes a two dimensional average unit cell. Likewise, the structure factor for a linear 
combination nk Q + mqo, n, m e Z is given by the (n, m) Fourier mode of P ao ^ (u, v). This 
means that the average unit cell, calculated for the wave vectors of the main structure and 
its modulation, can be used to calculate the peak intensities of any of the main reflections 
and its satellites of arbitrary order. Using (|T2j) and its generalization, it is possible to 
calculate the intensities of all peaks observed in the diffraction patterns. 



2.2 Patterson analysis 

Phase inversion problem pushed to development of a variation of the above method called 
average Patterson function analysis. 

It is well known that intensity I(k) can be obtained as a direct Fourier transform of 
the autocorrelation function 

g{x)= f°° P (e)p(e + *)d£ = /' x (YfM- >-,)){yh^ ■ .r -.r,))<K 

J— OO J— OO • 7. 



^2fjfkS(x- (x k -Xj)). 

3, k 



(13) 
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Let's suppose for a while that we deal with only one type of the atom (fj = f =const) 
in the structure. To calculate a Fourier transform of G(x) for scattering vector k we can 
utilize average unit cell formalism as well: again by introduction of the reference lattice with 
periodicity A = 2ir/k. Average Patterson function of the structure Sa(^) is a distribution 
of values of all interatomic distances reduced to average unit cell. We may however first 
convey x^ and Xj separately to the unit cell and subsequently compute the distribution - we 
may always add/subtract the multiplicity of the reference lattice to the distance; therefore 
Q\(u) may be regarded as the distribution of distances between these points' images in the 
cell. Since displacements sequences of Xk and Xj can be treated as independent variables 
(the last sum in (jlHj) is over all possible pairs of atoms) Sa is essentially an autocorrelation 
of probability distribution P\(u): 

/oo 
P x (Z)P x (Z-x)d£. (14) 
-oo 

Sa has all properties of P\; in particular we can use unreduced average Patterson functions 
where applicable and construct many-parameters average Patterson unit cells. 
The intensity: 

/oo 
Q\(x) exp(ikx)dx. (15) 
-oo 

2.3 Reference sites 

Let us now consider decorated structures. We will adopt two very reasonable assumptions 
about atomic structure which can be observed in nature. First, every relevant physical 
system consists of finite number of types of atoms and we will denote this number by 
N. Second, the atoms tend to organize themselves on scales comparable to interatomic 
distances while in the large scales every atom should be seen as distributed uniformly, close 
to its average concentration (we do not consider such effects like precipitations). Note 
that these assumptions allow us to work (for example) with every decorated quasicrystal 
obtained from projection method and also with classical crystals. 

We introduce idea of decoration sites, x n G A , where A fulfills assumption made upon 
A in the previous section. Each decoration site has a given type, denoted by the first 
capital letters of alphabet (A, B, . . . ). Relative numbers of given types of sites are (px, 
where X stands for A, B,. . . (p's are normalized: 

5>x = l- (16) 

x 

To every type of site we ascribe a group of atoms, with well defined scattering powers and 
positions in relation to the site; number of atoms in such group is rix (again X stands for 
A, B,. . . ) and 



(17) 



what means that atoms with the same scattering powers are considered to be different type 
if they are ascribed to sites with a different type. We will refer to the given atoms using 
notation Xj (J — 1, . . . , nx), what reads: j th atom in the group corresponding to site X. 
Scattering powers of atoms are f Xj and their relative positions with respect to the site a Xj • 
Classical crystallography uses only one type of sites organized in one of the Bravais 
lattice - the above idea is a generalization of classical picture. 

2.4 Overall displacements density function 

Displacement density function for X-type sites is CP^"(it). Possibly extended DDF for the 
atom of type Xj is simply 

P? i (u) = y*(u-a Xj ), (18) 

what follows from Theorem^ where we put /3 n — ax = const. We will refer sometimes to 
DDFs as "images" . As a result we get that the cell image of the given type of the decorating 
atom is just shifted image of the given type of reference site. The main advantage of working 
exclusively in the physical space can be easily^ seen. 

Now we are able build overall (E)DDFs P\(u), which Fourier transform given by (@J) 

will give us the structure factor for the whole decorated structure. In order to achieve that 

~x ■ 

we must add all P x 3 (u) weighted with scattering powers and relative number of given type 
of atoms: 

P ^= (EE^^a(«-«x,))M (19) 

X j=l 

where the normalization is 

n x 

^ = EE^^- ( 2 °) 

x j=l 

Reduced P\{u) constitutes average unit cell for the whole decorated structure and 
structure factor for scattering vector k = 2irm/X is then: 

/oo 
P x (u) exp(iku)du. (21) 
■oo 

2.5 Decorations in the language of Patterson function 

Usage of average Patterson function is crucial during the structure refining based on ex- 
perimental data due to the fundamental problem of phase lost. 
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According to (fUjl Sx in an autocorrelation of Pa- Written explicitly: 

^) = ^EEEE/^™ ^(o^e-^, (22) 

X Y j=l k=X J ~°° 

where X and Y run independently through A, B,. . .We can see that in the average cell 
N 2 components P x 3 * P x k may be pointed, each corresponding to a convolution of a given 
pair of decorating atoms' images. Intensity is given by: 

/oo 
5x(x)exp(ikx)dx. (23) 
-oo 

If I(2nm/X), n G Z, are known, we can rebuild Sa(^) by means of Fourier series as 
follows: 



5x(x) = y^Yl H^m/X) expH27rmx/A). (24) 

meZ 

Generalization to the case with many-parameters AUC is straightforward (compare to 
()12))): we give it here for 2 parameter unit cell: 

Sa, b (u,v) = T^r^ ^ ^2 I( y 2ir( y n 1 /a + n 2 /b))exp(-i2iT(n 1 u/a + n 2 v/b)). (25) 



3 Example: ID modulated crystals 



In this section we are going to present the way how the formalism works for the class 
of ID incommensurately modulated crystals. The method is of course not limited to ID 
structures and the example was chosen for its clarity. Generalization to real (3D) structures 
is possible (as suggested in the very beginning): we always deal with effective ID structures 
since structure factor contains scalar product of scattering vector and position vector, what 
corresponds to projection of the structure onto a given direction. Preliminary work shows 
that calculations for real 3D structures are characterized only by more complicated shapes 
of the distributions. 

Let us consider an aperiodic Sturmian sequence of two (positive) distances A and 33 
which sequence is coded using population vector p m = (p^iPm)' Pm+Pm = m an d Pm 
stands for number of X-type elements in the first m terms of the sequence: 



fm 



2 
Pm 



m 



a + 1 



m 



a 



(26) 
(27) 



where a e (0, 1) n (R \ 1 
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3.1 Decoration sites 

Set of decoration sites Ao contains points of the form: 

x m = p m -(A,-B). (28) 

We will say that the decoration site is X-type if x m+ i — x m = X, i.e. if the next site is 
in the distance X. 
Since 

A 

lim % = 1/a, (29) 

fm 

the relative concentrations of the given type of sites are 

^ = rb' ^ = rb- (30) 

One can easily recognize that A is just a model set obtained by projection a part of 
two dimensional square lattice onto a line with a slope tan _1 (cr). 
We will take our basic distances to be equal to respectively 

A = K ^ + a K , 'B=A/k, (31) 
k + a 

k G 1R_|_\{1}. a determines the sequence of A and 3 and k their relative length. 
It can be shown that reference sites' positions are given by: 

x m = ma + UM bQ (ma ), (32) 

where: 

M bo (ma ) = i - {x/b + (33) 

b =(a + l)a , U = A-'B, (34) 

and {x} = x — \\x — ||| is the fractional part of x. 



3.2 Average unit cells for decoration sites 

Form of eq. (}32|) suggests, that our structure is incomensurately modulated crystal with 
main period ao and period of modulation bo. For such a structure it has been shown (e.g. 
|0j), that diffraction pattern is supported on a discrete set of fc's of a form: 

k = n l k + n 2 q 0j (35) 

where ko = 2%/ao and go — 27r/6o- To calculate the whole pattern we will have to construct 
2-parameter average unit cell 7 a0j b (u,v), based on these scattering vectors. 
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Let's first build one parameters distributions. We can obtain that possibly unreduced 
displacement sequences of the decorating sites positions for k and q are respectively (we 
utilize Theorem P): 



a 

•Em 



bo 



UM bo (ma Q ) 



where : 



e = 1 - 6o/t7 = (1 + «r)/(l - «). 



(36) 
(37) 

(38) 



Using Kronecker's theorem j2j we can state that both of these sequences are uniformely 
distributed and hence corresponding (E)DDF's have just rectangular shape. 

It is also easy to show, that given types of sites reduced to average unit cell always 
occupy its well defined parts, as shown on Figure d Parameters of theses distribution are 
U and V = for k and qo respectively. 



U A = U- 



u B = u- 



a 



1 + a 

V A = ^— - — , V B = V 
1 + a 



1 + a 
a 

1 + a 



(39) 
(40) 



a) 
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b) 
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> 



P(u) 



A - type 
B — type 



Figure 1: Average unit cells for decoration sites, a) & b) show (possibly unreduced) AUC's 
for ko and go respectively. Given types of sites always occupy well defined parts of the cells. 
Picture shows situation for £ < 0. The distributions are centered around 0. 
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3.3 Overall density function 



Now let's consider 2-parameter average unit cell. Strong correlation between displacements 
sequences u' m and v' m given by (pT7jl makes the 2-parameter distributions to be supported 
on fragments of straight lines only. Each decorating atom brings additional component to 
the overall function of the form (compare to (fTHlH : 

P S,i>o( u > v ) = ^ 6 ( v - & - a ^ - Ux, c x + a Xl ). (41) 

U u U 1 , An . 

CA = -2TT^> CB= 2TT^ (42) 

and R(x;w,Xq) = 6(x — x + w/2) — 9(x — xo — w/2), where 9(x) stands for Heaviside 
step function. R is rectangle with center xq, width w and heigth 1. Example for 1 atom 
decorating A-type sites is shown in the Figure El 




Figure 2: Two parameter average unit cell for structure with 1 A-type atom decoration; 
also the positions of decoration sites were added (segments DE and EH are occupied by A 
and .B-type decoration sites respectively). Each decorating atom's image is an additional 
component fx l y?xPa l b { u i v )- I n the discussed case segment LK is an image of the dec- 
orating atom. CC = (dApCiAi) ~ shift vector. C = (ca,£ca), D = (—U/2,—£ i U/2) and 
F = (0, 0^(1 — £))• Such behavior is a consequence of image's shift in ID unit cell (inner 
panel), please compare to equation (|THj) . 

Once we have figured out geometry of the cell overall two parameter displacement 
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density function can be written using (fT§|): 



x x 

PooJboM = (EE^^4»(^))M ( 43 ) 



X 1=1 



where normalization 91 is given by (|20|). 

Structure factor for fc = nifco + ^2<?o is given by: 



" X 



F(k) = Y J Y,f^xF Xl} (44) 



x 1=1 



where: 



00 poo 

X, 



Fxi = I I P ao [ bo (u,v)exp(i(n 1 k u + n 2 q v))dvdu 

jka Xl ( exp (^(c x + U x /2)) - exp(iK(c x - U x /2))) 



oo J — CO 

i 



KU 



x 



sinf w ) 

= exp(i(Kc x + ka x ,)) — — ■ (45) 

w 

K = nik + ^n 2 go = k — n 2 q s , q s = qo(l — £) and w = KUx/2- Equation (jUjl is valid 
only for values of k given by f)35|) and allows us to calculate intensity of each peak in the 
pattern. We can treat (jUj) also as continuous function of scattering vector k, parametrized 
by the order of the satellite n 2 . In the latter case we obtain so called envelope functions 
E m (k), going through all satellites of the order m = n 2 . 

3.4 Comment on factor e lkax i 

When we add an atom at position xq with respect to the decoration site, the corresponding 
part of DDF in 2-parameter AUC is will be just the shifted image of this decoration site; the 
shift vector reads (xq, Xq). Partial structure factor of this shifted DDF for k = ni&o + n 2 q 
is: 

/OO /"OO 
/ P{u — x , v — x ) exp(i(nik u + n 2 q v))dvdu 
•00 J— 00 

P(u', v') exj>(i(niko(u' + x ) + n 2 q (v' + x ))dv'du' 

—00 j —00 

OO POO 

jkx I / p( u ^ v ') exp (i( ni k oU ' + n 2 q v'))dv'du'. (46) 



OO J —CO 
CO POO 



— CO J —CO 



The last integral is of course proportional to structure factor for not shifted DDF. The 
relevance of the proportionality factor will be discussed more deeply in the Summary. 
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3.5 Patterson function 

As we have mentioned the Patterson function is of the fundamental importance. General- 
izing equation (J22J) to 2-parameter AUC we obtain: 

~ 1 n x n Y _ _ 

5ao,b (u, V) = — ^ Yl Yl fXjfY k <Px<PY(P£ tbo * P^J, (47) 
X Y j=l k=l 

where the convolution P„ \ * P„ k , is given by: 

j poo poo 

— — / / 8{p-&-a Xl (l-S))5(p-v-S(a-u)-a Yk (l-ti))x 

UXVY J-oo J-oo 

R(a; U x ,c x + a x )R(a - u; U Y , c x + a x )dad/3 = 
=5(v - £u - (1 - 0{a Xl - a Yk ))T{u- 1, w, h). (48) 

T(u; I, w, h) = R(a; U x , c x +a X[ )R(a—u; U Y , c x + a Xl )da is a trapezium-like function; 
in our case I = \ \U X \ — \U Y \\, w = \U X \ + \U Y \, h = 2mm(\U x \, \U Y \)/(l+w). T was depicted 
in the Figure El 




Figure 3: Trapezium-like function used in equation (|48p. Average Patterson unit cell 
contains of sum of such functions with different parameters supported on segments of 
straight lines. 

Average Patterson unit cell contains N 2 components of the form given by (|48j): as 
ordinary average unit cell Patterson function function is supported on the fragments of 
straight lines only with slope £, but has more complicated shape. 

3.6 Pattern analysis 

In this section we present the algorithm allowing recovery of the structure if the diffraction 
pattern J(27r(rii/ao + ^2/^0)) i s known and we assume the crystal has internal structure as 
discussed above. 
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The pattern itself gives us values of k and q , and hence a and b ; using (|3*3)l we 
can compute a. Using we can inverse the transformation and obtain 9a ,&o( M ' 
Immediately a value of £ can be found, and subsequently using relation (J38|) the ratio 
k = A/H. Subsequently formulas (|3*T|) give values of A and 23. 

Average Patterson function can now be used to determine number of atoms, their types 
(scattering powers) and positions (type of decoration site and relative positions); of course 
due to phase problem we cannot avoid certain ambiguity. Detailed (numerical) algorithm 
will be discussed elsewhere. Many information can be extracted form average Patterson 
function relatively easily: note for example that from f)48j) . e.g. by investigating points 
where Patterson function crosses axes u or v, we can obtain information about relative 
positions of decorating atoms (ax k — ayj. 

Since average Patterson function is usually known with serious error, further refinement 
of atomic positions can be made by utilizing the structure factor. We can try direct fitting 
of the decorating distances including all of the diffraction peaks at once (e.g. by means of 
least squares method). It is also possible to use only subset of the pattern and optimize 
the fitness of envelopes E m {k). 

As it was mentioned earlier structure factor can be described as a sum of components 
F Xl given by which are products of two parts. The first one is the discussed phase 
factor, which depends only on relative decorating atom position; the second one is char- 
acteristic for the structure and does not depend on decoration. It is important feature 
which may improve the numerical algorithm. The part independent on decoration has to 
be calculated only once and it is enough to limit the refinement to the phase factors. 

3.7 Specific example 

As an example we present a structure with the following parameters: a = it, k = \[2 
and values of A and "B are chosen such that ao = 1 + 1/t 2 , r is a golden ratio equals 
(1 + y/E)/2. The reason for such strange-looking choice is to make correspondence with 
other publications; conventionally for er = 1/r and k = a we get famous Fibonacci chain. 
Above choice determines the value of other parameters, gathered in Table 13.71 

Decoration consists of two identical atoms with scattering powers equal to 1: the first 
occupies the the A-type segment in position = 0, the second B-type segment in position 
a Bl = S/2 « 0.628. 

Average 2-parameter Patterson function was plotted in the Figure 0] Please refer to 
the caption for detailed description. Figure presents diffraction pattern of the obtained 
structure with envelope functions. 

4 Summary 

The statistical approach to quasicrystals' decoration which has been presented in this paper 
can be intuitively understood as extension of the classical terms to the realm of non-periodic 
orderings. Average unit cell is extension of the classical unit cell; decorating atoms are 
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Value 



parameter 



exact 



approximate 



A 
B 

bo 

A: 
5« 




1.777 
1.256 
5.724 

-9.999 
4.547 
1.098 

12.074 



A/k 
(1 + tt)(1 + 1/t 2 ) 



1+7T 



2tt/(1 + 1/r 2 ) 

V(i + <7) 

90(1-0 



Table 1: The parameters of the system, which was used as the example. 



represented not as 5-like objects in the cell, but correspond to non-discrete "images" . Once 
we allow such generalization of standard notions, we can use the new concepts almost in 
the same manner as we used to work with the unit cell in the strictly periodic case. 

Structure factor is constructed in the way similar to classical case, as a sum of Fourier 
transform of similar objects with adequate phase factors (exp(ikax t ), for the example 
discussed), corresponding to atoms' positions. Such explicit separation of the decoration 
dependence from general properties of the structure can be a help during the refinement 
process. 

Simultaneous usage of average Patterson function and structure factor during analysis 
of the patterns can split the process into two separate parts. Patterson function gives us 
information about structure properties independent of the decoration. It can also provide 
us with number of decorating atoms and types of sites they are adherent to. Positions 
refinement can be done using structure factor (aided by the envelope functions). 
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Figure 4: Average Patterson function for the discussed example. Panel a) shows support 
of the function: it is a set of line segments with common slope £. Support crosses u axis 
for u = and u ~ ±0.691, what corresponds to all possible values (1 — — — 

±|ao (in this case). Panel b) presents perspective view of the cell itself. Results are 
numerical (histogram for about 20000 atoms in the structure) and fully supports analytical 
calculations. 
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Figure 5: Diffraction pattern for presented example. The envelope functions of orders 
m — 0, 1 and 2 were plotted. Pattern comes from numerical experiment, while envelopes 
were calculated analytically; full agreement exists for all diffraction peaks with given m. 
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